Spin-dipole strength functions of 4 He with realistic nuclear forces 
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Both isoscalar and isovector spin-dipole excitations of 4 He are studied using realistic nuclear 
forces in the complex scaling method. The ground state of 4 He and discretized continuum states 
with J n — 0~, 1~, 2~ for A = 4 nuclei are described in explicitly correlated Gaussians reinforced 
with global vectors for angular motion. Two- and three-body decay channels are specifically treated 
to take into account final state interactions. The observed resonance energies and widths of the 
negative-parity levels are all in fair agreement with those calculated from both the spin-dipole and 
electric-dipole strength functions as well as the energy eigenvalues of the complex scaled Hamiltonian. 
Spin-dipole sum rules, both non energy-weighted and energy-weighted, are discussed in relation to 
tensor correlations in the ground state of 4 He. 
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I. INTRODUCTION 

Spin-dipole (SD) excitations of nuclei have attracted 
much attention because of their connection with, for 
example, tensor correlation, neutron skin-thickness and 
neutrino-nucleus scattering. Especially the neutrino re- 
action involving light nuclei is important to the nucle- 
osynthesis at various stages. In the final stage of a core 
collapse supernova, the nuclei are exposed to the intense 
flux of neutrinos, and the neutrino-nucleus reaction rate 
is determined by the nuclear responses to such operators 
of the weak interaction as Gamow- Teller (GT), dipolc, 
SD, and so on [H, The SD operator brings about the 
first-forbidden transition of the weak interaction. In the 
case of N = Z nuclei, the allowed transition probabil- 
ity due to the weak interaction is small, and thus the 
first-forbidden transition can be a leading order, making 
a primary contribution to the cross section. 

Though the neutrino-nucleus reaction cross section can 
not be measured to good accuracy in a laboratory be- 
cause of its too small reaction rate, information on the 
spin excitation of nuclei can be obtained using a charge- 
exchange reaction. For example, in the (p, ri) or (d, 2 Hc) 
reaction at intermediate energy, the cross section at 
degree is a useful probe to extract the GT strength [3[ 
as well as the SD strength 0, @. The SD excitation 
can be obtained by measuring the cross section at larger 
angles 0, 0] ■ For a doubly closed shell nucleus the SD 
contribution is fairly large even at degree because of 
the hindrance of the GT strength. Much effort has been 
devoted to measure the SD transitions in light nuclei. 
Recently the charge- exchange reaction of ( 7 Li, 7 Be7) is 
undertaken to measure the electric dipolc (El) and the 
SD resonances of 4 He and 6 ' 7 Li [1, Q. A more recent 
measurement of polarization transfer observables with 
16 0(p,n) 16 F reaction [l(J indicates that valuable infor- 
mation on the SD excitations is attainable. 

The SD transition is also interesting in comparison 
with the El transition. The SD operator can change 
the spin wave function of the ground state, whereas the 



El operator can not. Since the SD operator has three 
possible multipoles, the study of its transition strength is 
expected to be more advantageous to see the spin struc- 
ture of nuclei than that by the El operator That 
is, this multipole dependence of the SD operator may be 
used to probe the role of non-central forces, especially 
the tensor force. The effect of the tensor force has in fact 
been studied theoretically by looking at the SD excita- 
tions in the shell model [12j and in the random phase ap- 
proximation based on Skyrme-Hartree-Fock [l3l [l4| and 
relativistic Hartree-Fock methods [l5j]. All these calcu- 
lations employed the variety of effective interactions and 
found that the residual tensor terms added to the inter- 
action play some multipole-dependent effects on the SD 
strength functions. 

The purpose of this paper is to study the SD excita- 
tions of 4 He using realistic nucleon-nucleon interactions. 
Only the ground state of 4 He is bound among A = 4 nu- 
cleon systems and its basic property is now understood 
fairly well thanks to several accurate methods for solving 
bound state problems of few-body systems All the 
excited states of 4 He are in continuum and its negative- 
parity states below the 2n + 2p threshold have J = 0, 1, 2 
with T = and 1. These resonances as well as the con- 
tinuum states may be reached by the SD operators. It is 
therefore quite challenging to accurately predict the SD 
strength function as a function of excitation energy be- 
cause we have to deal with the continuum states where 
not only two but also three particles may play an impor- 
tant contribution. On top of that we have to take into ac- 
count both short-range and tensor correlations due to the 
realistic nuclear forces [I?], [Hi • Very recently the present 
authors and Arai have done an ab initio calculation for 
the photoabsorption of 4 He [l9| using square integrable 
(C 2 ) basis functions in the framework of the complex 
scaling method (CSM) and have reproduced most of ex- 
perimental photoabsorption cross section data up to the 
pion threshold. A theoretical approach employed in the 
present paper is similar to that of the El case. 

It is well known that the ground state of 4 He contains 



2 



the £)-state (or the total spin S = 2 state) probability 
by about 14%, which is of course due to the tensor force. 
As shown in the calculation of bound-state approxima- 
tion [2^, [2l[ , the tensor force plays a vital role in correctly 
reproducing the spectrum of the excited states of 4 He. If 
one uses such effective interactions that contain no ten- 
sor components, there is no way to account for the level 
splittings of the negative-parity states of 4 He. Therefore 
the use of realistic nuclear forces is absolutely necessary 
for studying the SD strength in 4 He. In the same context 
we also study the charge-exchange SD transitions from 
4 He, leading to the negative-parity states of 4 H or 4 Li. 
We will pay due attention to the effect of the tensor cor- 
relation on the SD excitations. It should be noted that 
the SD excitation is here described based on the accu- 
rate ground-state wave function of 4 He @, [2l[. We also 
note that this study will serve fundamental data on the 
neutrino- 4 He reaction cross section in stars by integrat- 
ing the SD strength functions weighted by the neutrino 
energy distribution produced by the core collapse star. 

In Sec. [IT] we present our method of evaluating the SD 
strength functions, the CSM (Sec. Ill Ap and the C 2 basis 
functions (Sec. Ill Bp that are keys in the present calcu- 
lation. We show calculated results on the SD strength 
functions in Sec. Mil The SD strengths calculated from 
continuum-discretized states are presented in Sec. IHI Al 
The SD strength functions of both isovector and isoscalar 
types are displayed in Sec. IIIIBI A comparison of the 
peaks of the SD strength functions with the resonance 
properties of 4 He is made in Sec. lIII CI The SD sum rules, 
both non energy-weighted and energy-weighted, are dis- 
cussed in Sec. IIII Dl Conclusions are drawn in Sec. IIVI A 
multipole decomposition of the SD non energy-weighted 
sum rule (NEWSR) is discussed in Appendix A, and a 
method of calculating its relevant matrix element with 
our basis functions is briefly explained in Appendix B. In 
Appendix C we derive a formula that makes it possible to 
calculate the contribution of the kinetic energy operator 
to the SD energy- weighted sum rule (EWSR). 



II. CALCULATION METHOD OF SPIN-DIPOLE 
STRENGTH FUNCTION 

A. Complex scaling method 

The SD operator with the multipolarity A and its pro- 
jection fj, is defined by 



with 



Pi = n- x N , 



(i) 



(2) 



where r, is ith nucleon coordinate, xjy is the center-of- 
mass coordinate of the A^-nucleon system, and er^ is zth 



nucleon spin. The center-of-mass motion is completely 
removed in the present paper and only the intrinsic exci- 
tation is considered. The square bracket [pi x er^^ de- 
notes the angular momentum coupling of the two vectors 
or more generally the tensor product of spherical tensors 
to that operator specified by A fi. The value of A can take 
0,1, and 2. The superscript p of O^^ or Tf distinguishes 
different types of isospin operators, isoscalar (IS), isovec- 
tor (IVO), and charge-exchange (IV+ and IV—), that is, 
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is 



1, Ti 



ivo 
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IV± 



t±(i) 



(3) 



In the inelastic neutrino-nucleus reaction, the neutral 
current induces the IVO type operator as well as the IS 
one. The isospin operator t + = t x + it y (t_ = t x — it y ) 
converts a proton (neutron) to a neutron (proton) , which 
corresponds to the charge-exchange process X(n,p)Y 
(X(p,n)Y). 

The strength function of an initial state ^>q for the SD 
operator is defined as 

S( P ,X,E) = (* f \OL |*„> \ 2 5{E f -E -E) 



= -ilm5> |< 



A ^ E + E Q - H + ie 



0£ M l*o>, (4) 



where Sf^ represents a summation over fj, and all the fi- 
nal states tyf. Both the initial and final states are the 
eigenfunctions of a Hamiltonian H with the energies Eq 
and Ef. They are normalized as usual: ('IV I 1 ?,,) = <5„/„ 
and 5(E v i — E v ) for bound and unbound states, respec- 
tively. In the second expression of Eq. Q the summa- 
tion over the final states with the energy conservation of 
S(Ef — Eq — E) is converted to the imaginary part of a 
resolvent R 



R 



1 



E + Eq — H - 



(5) 



In the present paper we use the CSM to obtain the 
strength function. The CSM is widely used not only 
in atomic and molecular physics [22I [23j but in nuclear 
physics [24| as well. Very recently it has successfully been 
applied to calculate the photoabsorption cross section of 
4 He with a realistic Hamiltonian (l9j . The CSM allows 
us to obtain the strength function using only C 2 basis 
functions exclusively, making it possible to avoid an ex- 
plicit construction of the continuum state. The key of the 
CSM is to rotate both the coordinate and the momentum 
by a scaling angle 9 



Pj e 



-ie 



(6) 



which makes the continuum state damp at large dis- 
tances within a certain range of 9. The strength function 
S(p,X,E) reduces to 

S(p,X,E) 

= --Im]T (* | 0HU-H9)R(9)U(9)0% |* ) , (7) 
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where U(8) is the scaling operator that makes the trans- 
formation ^ and R(9) is the complex scaled resolvent 



R{0) = U{9)RU- 1 {8) = 



1 



E + E - H{0) + it 
with the rotated Hamiltonian 

H(9) = U{6)HU- 1 (6). 



(8) 



(9) 



Provided the eigenfunctions of H(9) are made to damp 
at large distances, they can be expanded with a set of C 2 
basis functions $>i(x) 



* v {e)=Y J c v i m i {x). 



(10) 



The coefficients C" (9) together with the complex eigen- 
value E v {0) are determined by diagonalizing H(9): 



H{9)^ v (9) = E V (P)V V (0). 



ill) 



The strength function S (p, A, E) is then calculated from 
the following expression: 



S(p,\,E) = ~J2 lm 



gg)gg(g) 
;r<— E + E - E u {6) + ie 



(12) 



where 



with 



v^(e) = {(* v (9)y\oi fl (6)\um ), 

= < W)*o)*| |<M*)> (13) 



0^(0) =<V- ^,W=<e^. (14) 



Note that U(9)^o is here taken to be the solution of Eq. 
(fTTj) corresponding to the initial state. 

If a sharp resonance exists, the angle 9 has to be ro- 
tated to cover its resonance pole on the complex energy 
plane [23j, [24( . Practically the scaling angle 9 is chosen 
by examining the stability of the strength function with 
respect to 9. See Refs. [l^, [2l| for some examples on the 
9- dependence. 



B. Correlated Gaussians and global vectors 

1. Hamiltonian 

The Hamiltonian H we use contains two- and three- 
nucleon interactions 

N 

H = 2J Ti - T cm + ^2 Vij + X! Vi 3 k ■ ( 15 ) 

i— 1 i<j i<j<k 

In the kinetic energy Ti the proton-neutron mass differ- 
ence is ignored. Two different two-nucleon interactions, 



AV8' |25[ and G3RS [26| potentials, are employed to ex- 
amine the extent to which the strength function is sen- 
sitive to the -D-state probability of 4 He. The L 2 and 
(L ■ S) 2 terms in the G3RS potential are ignored. The 
AV8' potential is more repulsive at short distances and 
has a stronger tensor component than the G3RS poten- 
tial. As the three-body interaction (3NF) we adopt the 
spin-isospin independent phenomenological potential 27] 
that is adjusted to reproduce both the inelastic elec- 
tron scattering form factor to the first excited state of 
4 He as well as the binding energies of 3 ' 4 He and 3 H. 
The Coulomb potential is included, but the isospin is 
treated as a conserved quantum number. The nucleon 
mass mjv and the charge constant e used in what follows 
are h 2 /m N = 41.47106 MeVfm 2 and e 2 = 1.440 MeVfm. 



2. Basts functions for bound states 

We solve the four-body Schrodinger equation using a 
variational method. A choice of the variational trial func- 
tions is essential to determine the accuracy of the calcu- 
lation. A bound-state solution with spin-parity J n of 
iV-nucleon system may be expressed in terms of a linear 
combination of the LS coupled basis functions 

^LS)JMjTUt = A [ ( t ) lx Xs] JMj VTMt , (16) 

where A is the antisymmetrizer, and the spin function 
Xs is given in a successive coupling as 



XSl2,Sl23,---,SMs 

= [-[[Xi(l)xn(2)k 2 xn(3)] Sl2 



\SM S 



(17) 



Note that the above spin function forms a complete set 
provided all possible intermediate spins (S12, £123 j ••■) 
are included for a given S. The isospin function t\tm t is 
given in exactly the same way as the spin function. 

The spatial part <f>\ should be flexible enough to cope 
with the strong tensor force and short-range repulsion. 
The tensor force mixes the S and D components in the 
wave function and the short-range repulsion makes the 
amplitude of the two-nucleon relative motion function 
vanishingly small at short distances. Many examples 
show that the correlated Gaussian (CG) bas is 1281. 1291 is 
flexible enough to meet these requirements [lg, [20|, \3C_ 
See a recent review 31] for various powerful applications 
of the CG. Let an (N — l)-dimensional column vector or 
an (N — 1) x 1 matrix x denote a set of relative coor- 
dinates whose «th element is a 3-dimensional vector 2^. 
A set of the Jacobi coordinates is most often employed 
for x but other sets of relative coordinates may be used 
as well. The spatial part given in the CG with two 
global vectors (GV), takes a form (32l - [35j 



F( Li l 2 )lm l (ui,u 2 ,A,x) 
= exp{-\xAx) [yL l {uix) x yh 2 {u2x) 



LMl 



(18) 
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with 



yim{v) = v e Y £m (v), 



(19) 



where A is an (N — 1) x (N — 1) positive-definite, sym- 
metric matrix and xAx is a short-hand notation for 
SfjUi AijXi ■ Xj. The tilde stands for the transpose of a 
matrix. Parameters U\ and u 2 are (TV — l)-dimensional 
column vectors that define the GVs, u\x{= X^ili. 1 unXi) 
and U2X, and these characterize the angular motion of the 
system. 

The CG-GV basis (fT8j) apparently describes correlated 
motion among the particles through the off-diagonal el- 
ements of A and the rotational motion of the system is 
conveniently described by different sets of (Li^Li) car- 
ried by the two GVs. Most noticeable among several 
advantages of the CG-GV basis functions are that the 
functional form of Eq. (|18[) remains unchanged under an 
arbitrary linear transformation of the coordinate x, that 
the matrix elements for most operators can be evaluated 
analytically, and that the formulation can readily be ex- 
tended to systems with larger N. Useful formulas for 
evaluating matrix elements with the CG-GV basis are 
collected in Appendices of Refs. [34l. l35|. 

All possible L, S sets are adopted to specify the ba- 
sis functions for a given J. The value of S can be 
0, 1, and 2 for the four-nucleon system, and all pos- 
sible L values that make J with S are included. For 
a given L* we choose the simplest combination of (£1, 
L2): (L\ — L, L2 — 0) for a natural parity state with 
7r = ( — 1) L and (Li — L, Li = 1) for an unnatural par- 
ity state with it = (— l) i+1 , respectively. An exception 
is that no basis function with L n = 0~ is included in 
our calculation because that special case needs at least 
three GVs [35|,[36j]. It should be noted, however, that the 
= CP configuration may be excited by the El and SD 
transitions only through the (L = S = 1) component of 
the ground state of 4 He. Since the probability of finding 
that component is quite small (less than 0.4%) [lH [34j. 
practically we do not miss any SD strength by the neglect 
of the L n = CP configuration. This is really the case in 
the El strength function [l9[ and in the SD case as well 
as shown in Sec. IIII Dl 

The parameters, A, u±, and 112, are determined by the 
stochastic variational method (SVM) [32|, HH, H3 ■ The 
calculated properties of 3 H, 3 He, and 4 He agree with ex- 
perimental three- and four-nucleon data very well [2l[ . 



Square-integrable basis functions for spin-dipole 
excitations 



We construct the basis functions for the final states 
with J~T that are excited by the SD operator of 
A = J. The accuracy of the CSM calculation depends on 
how fully the basis functions are prepared. In Ref. (l9j . 
the present authors and Arai described a way to con- 
struct the four-body continuum-discretized states with 



jwrp _ rp^e guidelines of the construction were to 

take into account both sum rule and final state interac- 
tions between the particles in the continuum. The total 
photoabsorption cross section is calculated via the El 
strength function and it succeeds to reproduce the mea- 
sured cross section up to the pion threshold. Here we 
take the same route as that of Ref. [19( with a possible 
modification due to the spin flip of the SD operator. 

We define a single-particle (sp) basis, which describes a 
single-particle like excitation from the correlated ground 
state by the SD operator. This class of basis functions is 
expected to play a vital role in accounting for all the SD 
strength. The basis is constructed as 



*7 



A 



x Xs> 



T ,xyx(pi) 



A/7 



VTMt 



(20) 



where </>J. (i) is the space part of the ith basis function of 
a truncated ground-state wave function, ^Jqoo' °^ 4 ^ e - 
The wave function ^q 000 consists of [<Pq x xo] 00^00 and 
[4>t x X2]oot?oo, with any configurations of [4>f x Xi]oo%o 
being omitted, which leads to 1.53 MeV loss for the 
ground-state energy of 4 He. See Ref. [19j for the detail. 
As for the spin part, differently from the El case [l9| we 
take into account the complete set for a given S' , which, 
depending on the total spin Si of the ith basis function 
of ^oooO' ^ s cnosen as S" = 1 for Si — and S' = 1 and 
2 for Si = 2, respectively. 

The 3iV + N two-body and d + p + n three-body disin- 
tegration channels arc defined in the same manner as in 
Ref. [l9|. 

The calculations are performed not only in each basis 
set of sp, 3iV + N and d + p + n but also in the 'Full' 
basis that includes all of them. The number of basis 
functions in the Full model with the AV8'+3NF poten- 
tial is 2980, 6400, 6540, 4380, 8800, 9540 for J n T = 
Q-Q, 1-Q, 2-Q, 0~1, 2~1, respectively. 



III. RESULTS AND DISCUSSIONS 

A. Discretized spin-dipole strength 

First we show the continuum-discretized SD strength. 
For this purpose the Hamiltonian is diagonalized in the 
basis states that are defined in Sec. Ill B 31 This calcu- 
lation corresponds to the CSM solution with 6 = 0°. 
Figure [T] displays the reduced transition probability for 
the IV SD operator (p = IV0) 

B(p, J-T, v) = K*t~ MT (6 = 0°)PU*o )| 2 , (21) 

Mil 

as a function of the discretized energy E v (6 = 0°), where 
v is the label to distinguish the discretized energy. Here 
A is equal to J and T = 1 . 

The results of calculation are similar to the El case 
of Ref. [l!| . In the calculation with the sp configu- 
ration only, the strength is concentrated at one state 
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FIG. 1: (Color online) The discretized isovector SD reduced transition probabilities of IV0 type for 4 He as a function of 
excitation energy. The JT values of the excited states are 0~1, 1~1, and 2~1, respectively. The transition probabilities are 
displayed, arranged vertically for each JT case, depending on the configurations included in the calculation. See text for 
details. The AV8'+3NF interaction is used. 



that appears at about 25 McV for all the cases with 
J V T = 0~1, 1~1, 2~1. The state may correspond to 
the observed level at 23.33 MeV for J"T = 2"1 and 
25.28 MeV for J n T = 0~1, respectively For the 

J^T = 1~1 case, two levels with very broad widths are 
known at 23.64 and 25.95 MeV. Since the energy of the 
prominent SD transition strength is lower than that ob- 
tained for the El transition strength [19|, the 23.64 MeV 
level probably has SD character, whereas the 25.95 MeV 
level is excited by the El operator. 

Similarly to the El transition strength fl9| . two or 
three peaks are obtained with the 3N + N configuration 
and relatively small strength is spread broadly above 30 
MeV. The prominent peaks below 30 MeV shown in the 
3N + N calculation continue to remain in the Full ba- 
sis calculation, which again confirms the importance of 
the 3N + N configuration to describe the low-energy SD 
strengths as in the El strength. We also calculate the 
SD strength with the G3RS+3NF interaction. Both dis- 
tributions look similar, indicating the weak dependence 
of the SD strength on the realistic interactions employed. 

The so-called softening and hardening of the SD exci- 



tation is discussed in Refs. [13j, [14|, where the residual 
tensor force is turned on or off and the energy and the 
strength of the SD excitation are compared each other. 
We think that the conclusion drawn by such comparisons 
is not always true because switching off the important 
piece of the nucleon-nucleon interaction may cause a sig- 
nificant change in the continuum structure that can be 
reached by the SD operator. In fact, we can not turn off 
the tensor force. If the tensor force were switched off, 
the ground state of 4 He would not be bound and more- 
over the spectrum of the negative-parity states would be 
far from the observed one [2(| Hl| . As will be discussed 
later, we find no quenching of the SD strength but con- 
firm that our SD strength calculated with the realistic 
nuclear forces satisfies the NEWSR perfectly. 

Figure [5] displays the reduced transition probability for 
the IS SD operator (p — IS). The gross structure of the 
strength is similar to the IV SD case. In the sp configura- 
tion calculation, strongly concentrated peaks ap pear at 
the energies not far from the observed levels [3g|: 21.01, 
21.84, 24.25 MeV for J n T = 0~0, 2~0, 1~0, respectively. 
The importance of the 3N+N configurations is indicated 
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FIG. 2: (Color online) The same as Fig. [TJ but for the isoscalar SD transitions. The JT values of the excited states are 
0~0, 1~0, and 2~0, respectively. 



by those peaks that appear in the 37V + N configuration 
calculation and continue to exist in the Full calculation. 
The Full calculation predicts one prominent peak at 20.85 
MeV for J^T = 0~0, which may correspond to the 21.01 
MeV level with the small decay width of 0.84 MeV [H|. 
In the case of J^T = 2~0, the two prominent strengths 
are obtained at the energies close to each other, suggest- 
ing a relatively small decay width. The relationship be- 
tween the strength functions and resonance properties 
will be discussed in Sec. MI CI 

Three negative-parity states of 4 He with J = 0, 1, 2 
and T = are observed slightly above the 2n + 2p thresh- 
old of 28.3 MeV 38]. These states may be excited by the 
IS SD operator. In fact, a comparison of the IS and IV 
SD strengths obtained in the d + p + n configurations 
clearly suggests that more strength is found in the IS 
case around the excitation energy of 30 MeV. Therefore 
some discretized states at around 30 MeV shown in Fig. [2] 
may be precursors of those observed states. Since the 
three observed states almost entirely decay by emitting 
deuterons, it is likely that they have d + d structure with 
a P-wave relative motion. The P-wave relative motion is 
possible only when the channel spin of two ds is coupled 



to 1 [35j. Thus we have the possibility of J n = 0~, 1~, 
and 2~ in accordance with the observation. Our basis 
functions partially include the d + d type configurations, 
but an explicit inclusion of them may be desirable to dis- 
cuss this issue in more detail. 

It is interesting to recall the enhanced SD excitations 
of the first excited state JT = 0^0 that is interpreted 
as having a well developed 3N + N ( 3 H+p and 3 He+n) 
structure [2(| It is shown in Ref. that some 

negative-parity states can also be understood as par- 
ity inverted partners of the first excited state and the 
SD transition strengths from that state are quite en- 
hanced and mostly exhausted by only those negative- 
parity states. Figure|3]exhibits the SD reduced transition 
probabilities from the 0^0 state as a function of excita- 
tion energy. The transition probabilities of both IS and 
IV0 are very much enhanced, approximately 20-30 times 
larger than those from the ground state and each of the 
strengths is concentrated at the respective peak. The 
excitation energies of the peaks are 20.85, 21.37, 21.30 
MeV for J^T = 0^0,1-0,2-0 and 21.10, 21.32, 21.33 
MeV for J n T = _ 1, 1~1, 2 _ 1, respectively. The energy 
required for the 0^0 state to reach the peak position is 
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FIG. 3: (Color online) The discretized SD reduced transition probabilities of IS and IV0 types from the first excited + state 
of 4 He as a function of excitation energy. The calculation is done in the Full basis using the AV8'+3NF interaction. The 
calculated excitation energy of the + state is 20.33 MeV from the ground-state of 4 He [2l| . 



only 0.5 — 1.0 MeV. The neutrino reaction rate would be 
greatly enhanced if there were such a situation in which a 
plenty of the first excited states of 4 He existed in the core 
collapse star. The situation may, however, be unlikely as 
the life time of that state is short and its excitation en- 
ergy (20.21 MeV) is considerably high compared to the 
typical temperature of the collapsing star [39( . 



B. Spin-dipole strength functions 
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FIG. 4: (Color online) Isovector SD strength functions of IV0 
type and El strength function for 4 He as a function of exci- 
tation energy. 



In what follows we will present results obtained in the 
Full basis calculation with the AV8'+3NF potential using 
the scaling angle 6 = 17° unless otherwise mentioned. We 



count the excitation energy of the continuum state from 
the calculated ground-state energy of 4 He that is listed in 
Table 1 of Ref. 21]. Preliminary results on the GT and 
SD strength functions were reported in Refs. (2ll. Eoj. 

Figure H] plots the SD strength functions of IV0 type. 
For the sake of comparison, the El strength function 
is also plotted by choosing the El operator as Mi^ = 

Y^Li PiuhQ ~ T z(*))- A s seen m the figure, the three SD 
strength functions show narrower widths at their peaks 
than the El strength function. Moreover their peak po- 
sitions including the El case well correspond to the ob- 
served excitation energies of the four T = 1 negative- 
parity states of 4 He [38(. We will discuss this point in 
Sec. Imci 

Figure [5] displays the charge-exchange SD strength 
functions of IV± type as well as the charge- exchange El 
strength that is excited by the operator 



M 



iv± 
in 



N 
i=l 



J, 



IV± 



(22) 



Since the mass difference between protons and neutrons 
is ignored in the present calculation, we need to shift the 
calculated energies of 4 H or 4 Li by ±(m„ — m p ). This 
adjustment makes it possible to correctly reproduce the 
thresholds of 3 H+n for 4 H and 3 He+p for 4 Li, respec- 
tively. Similarly to the IV0 case, the excitation energies 
of the charge- exchange SD peaks correspond to the ob- 
served levels of 4 H and 4 Li, and their widths are narrow 
compared to the charge- exchange El strength function. 

We display in Fig. [B] the IS SD strength functions that 
reflect the JT = A~0 continuum states of 4 He. These IS 
SD strength functions, especially for the 0~ and 2~ cases, 
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FIG. 5: (Color online) Isovector charge-exchange SD strength 
functions of types IV+ (a) and IV— (b) for 4 He as a function 
of excitation energy. The excitation energy is counted from 
the ground state of 4 He. 



show much narrower distribution than the IV strength 
functions. These peak energies again appear to corre- 
spond to the observed T = negative-parity levels in 
4 He. A close comparison between Figs. |6] and H] indicates 
that the 0~ case is noteworthy compared to the 1~ and 
2~ cases in that the energy difference in the peak posi- 
tions of the same J~ becomes much larger. As discussed 
in detail in Refs. [2(J [HI, the reason for this is under- 
stood by analyzing the role played by the tensor force 
among others. In the previous subsection, we mention 
the three negative-parity states with T = that are ob- 
served slightly above the four-nucleon threshold and are 
expected to have d + d structure. Though no concen- 
trated strength suggesting such states is seen in Fig. [51 
the falloff of the IS SD strength around 28 - 30 MeV 
looks flatter than that of the IVO case especially in the 
J n = 1~ state. This indicates that some IS SD strength 
may exist in that energy region. To be more conclusive, 
however, a study including d + d configurations explicitly 
is desirable. 

To the best of our knowledge, there are no data that 
can directly be compared to the theoretical strength func- 
tions presented above. An only exception is the measure- 
ment of the charge- exchange reaction 4 He( 7 Li, 7 Be7) [8], 
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FIG. 6: (Color online) The same as Fig.[4]but for the isoscalar 
SD strength functions. 
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FIG. 7: (Color online) Summed isovector SD strength func- 
tions for the process from the ground state of 4 He to 4 H as 
a function of excitation energy. Relevant experimental data 
taken from Ref. || are plotted for reference. See text for the 
details. 



9], from which the spin-nonflip (AS = 0) and spin- flip 
(AS = 1) components are separated by measuring the 
0.43 MeV 7 ray of 7 Be in coincidence with the scattered 
7 Be. The former cross section is ascribed to the El tran- 
sition, while the latter to the SD transition. The shape 
of the deduced photoabsorption cross section fairly well 
agrees with other direct measurements using photons (see 
Fig. 9 of Ref. [l9(), but the absolute magnitude is not de- 
termined definitively. The SD spectra corresponding to 
the excitation of the 4 H continuum from the ground state 
of 4 He is extracted from the spin-flip cross section in a 
similar way. Figure [7] compares the SD strength functions 
of type IV+ with the 'experiment'. In this figure the the- 
oretical curve represents just a sum of the strength func- 
tions with A = 0,1,2, and the experimental distribution 
is normalized in such a way that both strength functions 
give the same strength when integrated in the energy re- 
gion from E—18 to 44 MeV where the experimental data 
are available. The comparison between theory and exper- 
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iment in Fig. [7] should thus be taken qualitative as several 
assumptions are made in the analysis of the experiment. 
The peak observed at 24 MeV agrees with the calculated 
one (see also Fig. 03a)) and certainly it corresponds to 
the J^T = 2~1 resonance of 4 H. We see some difference 
in the shape of the strength function. Two conceivable 
reasons for it include firstly that the spin-nonflip process 
can in fact contribute to the SD transition as the ground 
state of 4 He contains 5 = 2 components and secondly 
that some higher multipole effects may contribute to the 
cross section particularly at high energy [l(J. The first 
reason is easily understood if we consider the transition 
from (L = S = 2) to (L = 1, 2,3, S = 2). Further exper- 
imental information is needed to make a direct compari- 
son with the calculation. 



C. Resonance parameters 

As noted in Sec. lIIIBl all the SD and El strength func- 
tions exhibit some common feature: They all have one 
peak, though the width of the strength distribution de- 
pends on the multipolarity A and the isospin T. It looks 
quite reasonable to identify the peak as a resonance. The 
resonance energy may be identified with the energy where 
the peak is located. We also estimate the decay width 
of the resonance by the difference of two excitation ener- 
gies at which the strength becomes half of the maximum 
strength at the peak, which agrees with a correct width 
if the strength function shows the Lorentz distribution. 
Actually the distribution is not Lorentzian in general as 
we see below, but this crude estimate should be useful as 
a guide. Table |l] lists the resonance energies and widths 
of the negative-parity states of 4 He, 4 H, and 4 Li that are 
determined in this way. The agreement between theory 
and experiment is very satisfactory. The average devia- 
tion of the calculated resonance energies from experiment 
is less than 0.4 MeV for 4 He despite the fact that most of 
their widths are larger than 5 MeV. The estimated width 
is also in reasonable agreement with experiment. 

A four-nucleon scattering calculation that couples 
3 H+p, 3 He+n, and d+d channels as well as many pseudo 
states is performed in Ref. [35[ using the same Hamilto- 
nian as the present study. Though the calculated phase 
shifts for the J^T = 0~0 state show a clear resonance 
pattern at the energy consistent with the 0~0 level of 
4 He, the phase shifts of the 2~0 and 1~0 states do not 
rise high enough to enable one to extract the resonance 
parameter. A more sophisticated analysis is needed to 
reveal resonances using, for example, the time-delay ma- 
trix [4l], ■ In this context we may say that extracting 
the resonance parameter from the strength function is ro- 
bust and can be applied to any case where even no sharp 
resonance is expected. 

Since the resonance parameter obtained above is not 
directly determined from the complex eigenvalue E^IO) 
of the Hamiltonian, one may argue that the agreement is 
fortuitous. Of course it would be very hard to predict the 



resonance parameter correctly if a chosen operator is such 
that has only tiny strength to that resonance. It is there- 
fore interesting and important to examine the complex 
eigenvalues that constitute the basis of the strength func- 
tion. To this end we rewrite the strength function (fT2|) 
as 

S( P ,X,E)--^ {E _ £m)2 + k{im)2 , 

(23) 

where e v {6), 7^(6*), a u {9), and j3 v {9) are defined by 
E v {e) = e v (e) + E - l - lv {6), 
E 3 V 0?) = (<0 + i% X (0) ■ (24) 

The first term of the numerator of Eq. (|23|) gives the 
Lorentz distribution, while the second term contributes 
to the background distribution. In principle a resonance 
may be identified as such E v (6) that is stationary with 
respect to the variation of 9 [23j. Then the strength 
function (|23[) has a ^-independent peak around such a 
stationary energy e v . Resonance parameters of electron 
and positron complexes are in fact determined very well 
by examining the ^-trajectory of E v {6) [H, HJ]. This 
is possible because H{9) for the atomic case has simple 
structure, H(9) = Te^ 2W + Ve~' w , where T and V are the 
kinetic energy and the Coulomb potential energy. In the 
nuclear case, however, H{9) is by far complicated and a 
large-angle rotation of the nuclear potential may lead to 
a very long-ranged potential, which, together with inher- 
ent difficulties in solutions with the nuclear Hamiltonian, 
makes an accurate solution of Eq. (|II[) extremely hard. 
Therefore, we first look for such eigenvalues that deviate 
from the rotating-continuum line as possible candidates 
for a resonance and choose the one that is closest to the 
peak energy of the strength function. 

For the JT = 0~0 and 2~0 states, which have a rela- 
tively small decay width, we find only one candidate that 
may correspond to the observed resonance but other J^T 
states have two or three candidates below 4A^ threshold. 
However, no candidate is found for the 1^1 state that 
is excited by the El operator. The resonance energies 
and widths determined in this way are also listed in Ta- 
ble |TJ The resonance energy obtained from the complex 
energy eigenvalue is in excellent agreement with experi- 
ment, even better than that determined by the strength 
function. The width is also satisfactorily reproduced. 
Two approaches to determining the resonance parame- 
ters produce successful results, and they are powerful, 
robust, and complementary. 

Figure [5] compares with experiment the resonance en- 
ergies of the negative-parity states of 4 He, 4 H, and 4 Li 
that are determined from the complex energy eigenvalues 
and the strength functions. It is striking that the theory 
reproduces the experimental spectrum in correct order 
and moreover closely to the observed excitation energy. 



10 



TABLE I: Resonance energies Er and widths T, given in MeV, of negative-parity levels of ^4=4 nuclei. Calculated values are 
extracted from the complex eigenvalues E(6), and the SD and El strength functions S(E). Experimental data are taken from 
Ref. 
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FIG. 8: Negative-parity levels of A=4 nuclei. Excitation energies are referred to the ground-state of 4 He. Solid, dashed 
and dotted lines of calculation are based on the complex eigenvalues, the SD and El strength functions, and bound-state 
approximation [2lj . respectively. Experimental data are taken from Ref. Levels belonging to the same J n T are connected 
by thin dotted lines. 



The dotted line in the figure denotes the energy obtained 
with a kind of the real stabilization method [45[, that 
is by diagonalizing the Hamiltonian in the CG-GV basis 
functions [2l| . Here the SVM search is performed to op- 
timize the parameters of the basis functions by confining 
the four nucleons in some configuration space. It should 
be noted, however, that such calculation faces difficulty 
when dealing with a resonance with a very broad width 
such as the 1^~ level of 4 He, and therefore the resonance 
energy obtained in that calculation should be taken only 
approximate. 



D. Spin-dipole sum rules 

Sum rules are related to the energy moment of the 
strength functions in different order and can be expressed 
with the ground-state expectation values of appropriate 
operators from which we can obtain interesting informa- 
tion on the electroweak properties of nuclei [SSL l47l | . 

Throughout Sec. MIDI and Appendices A and C we 
denote the numbers of nucleons, neutrons, and protons 
by A, N, and Z, respectively. Accordingly the center- 
of-mass coordinate is instead of xjy. In the other 
sections N is used to denote the number of nucleons be- 
cause the symbol A is reserved to stand for the matrix 
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that appears in Eq. ([18 



Non energy-weighted sum rule 



Here we discuss the NEWSR for the SD operator 



m (p,A) = / S(p,X,E)dE. 
Jo 



(25) 



The use of the closure relation enables us to express 
the NEWSR to the expectation value of the operator 
Eu^Au^Au w ^h rcs pcct to the ground state It is 
convenient to express that operator as a scalar product 
of the space-space and spin-spin tensors 



2? K) o = E 



[ Pi x Pj ] K ■ [vi x <r 3 ] K )TfT^ (26) 



where the rank k can be 0, 1, and 2, and the symbol 
(T K ■ V K ) = (-1) k v / 2k + 1[T k x V K ] 00 denotes a scalar 
product of spherical tensors, T K and V K . As shown in 
Eq. (JXTJ of Appendix A, the NEWSR ^Bl is equivalently 
expressed, with use of U\ K of Eqs. (|A3|) and (|A4I) . as 



m (p,\) = ^ ZTv, 



\Q(«)o) 



(27) 



The expectation value, (Q^q) 



(*o|Qf K)0 l*o>, can be 



evaluated using the basis functions (fl"6l) . ([IT]) , and ([T8]). 
as explained in Appendix B. 

In order to check the extent to which the NEWSR 
is satisfied, we compare too (p, A) that is calculated sep- 
arately with Eq. ([35]) or with Eq. ([27]). Table [H] lists 
the calculated NEWSR for the SD strength functions. 
We also list the values of (Q P K \ ) in Table [ill] for the 
sake of discussions below. As seen in Table HU the two 
different ways of calculating the sum rules give virtu- 
ally the same result for both cases of AV8'+3NF and 
G3RS+3NF interactions, which is never trivial because 
we use the fully correlated ground-state wave function 
for 4 He. The perfect agreement confirms that the basis 
functions prepared for the description of the SD excita- 
tion are sufficient enough to account for all the strength 
in the continuum. The NEWSR calculated with Eq. ([27]) 
for the Minnesota (MN) potential [48| is also listed in 
Table [IT] A comparison of the central MN force case with 
the realistic potentials will be useful to know how much 
the sum rule is affected by the tensor force. 

Among the three expectation values of (Q^q) m 
Eq. (1271) . the n — term gives a dominant contribution 
to the NEWSR. See Table [nil This is obviously because 
the major component of the ground state of 4 He is S = 
and it has a non-vanishing expectation value only for 
Q P ^ - In this limiting case mo(p, A) is proportional to 
U\o- Therefore the A-dependence of the NEWSR turns 
out to be 1 : 3 : 5 for A = 0, 1, 2, independently of p. 



This rule is confirmed in the MN case of Table [TT] The 
deviation from this ratio is due to the contributions of 



(2)0 



other Q P r K \ Q terms, especially the k — 2 term. The Q 
term contributes to the NEWSR through the coupling 
matrix element between the S — and 5 = 2 compo- 
nents of the ground state of 4 He. Since the admixture 
of the S* = 2 component is primarily determined by the 
tensor force, the deviation reflects the tensor correlations 
in the ground state. Neglecting the minor contribution 
of Q P ^ , Eq. ([2"T]) suggests that mo{p, A) is very well ap- 
proximated by 



3 v-(0) 

m (p, 1) = m Q (p,0) 



1 



m (p, 2) 



2 «Q[ 1)0 )-3(Q^ 2)0 » 
3 

2 \^(2)o)' 

rno(p,0) + kHQ p 1)0 )-3{Q P 2)a)) 



m o (p,0) 

5 
3 



:m (p,0)--(Qf 2)0 ) 



(28) 



Thus the deviation of the ratio from 1 : 3 : 5 is simply 
controlled by — | (Q P 2 } )i which is very well satisfied in 
the examples of Table [IT] Since (Q p 2 ) ) ^ s negative for 
p = IS, the ratio further increases from 1:3:5, whereas 
it is positive for p = IVO and IV±, and the ratio approx- 
imately reduces to 1 : 2 : 4. 

As discussed above, (Q P K \ ) plays a central role to de- 
termine the NEWSR for the SD strength functions. In- 
verting Eq. ([27]) makes it possible to express (Q^q) as 
a sum, over the multipole A, of the NEWSR 



(2^ )0 ) = E C/ ^o(p,A), 

A=0 



(29) 



where U~ x is the inverse matrix of U as given in Eq. (|A5j) . 
If the NEWSR for all A are experimentally measured, the 
above equation indicates that (Q^ k j ) for all k can be 
determined from experiment. Some examples are 

3 (2(o)o) = m °(P> °) + m o(P' 4 ) + m o{p, 2), 

6 (2( , 2)0 ) = 10m (p, 0) - 5m (p, 1) + m (p, 2). (30) 

To clarify the physical meaning of the operator Q P K ^ , 
it is instructive to decompose it into one- and two-body 
terms: 



qv _ qpW , qpC- 



-( re )0> 



(31) 



where 



eS = E (te x pA« ■ x 'A*) 7 * ( 32 ) 

j>i=i 



A 
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TABLE II: Non energy- weighted sum rules of the SD strength functions, in units of fm 2 , calculated from different models for 
the nucleon-nucleon potentials. The sum rules calculated by Eqs. (|25[) and (I27[) are labeled mo(p, A) and SR, respectively. 
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TABLE III: Expectation values of and its one- and two-body terms with respect to the ground state of 4 He. Values are 

given in units of fm 2 . 
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with 



rpp\rpp 

3 * ' 



(33) 



The isospin operators in Eq. (1321) are simplified with use 
of Eq. (A2): T[ f T[ is 1 for p = IS, IVO, and (l=Fr z (i))/2 
for p = IV±, whereas T?- is 2 for p — IS, 2r z (i)T z (j) for 
p = IVO, and ((t(<) • r(j)) - T z (i)r z (j))/2 for p = IV±, 
respectively. The one-body term is spin-independent and 
appears only for k — 0, which gives the largest contribu- 
tion to the NEWSR. The two-body term with k = 2 is 
particularly interesting because it contains the tensor op- 
erator characteristic of the one-pion-exchange potential. 
See Appendix A for detail. 

The expectation value of the one-body term is ex- 
pressed in terms of the root-mean-square radius of nu- 
cleon distribution in the ground state 



IV0(1) 
(0)0 

2 ) 

pi ' 



= A(r 

-)IV-(l) 



(0)0 



)=N(r 2 n ). 



(34) 



Noting that the two-body term Q™^ 2 * 1 is identical to 

Q(^)q for any k, we obtain the following well-known 
relation between the NEWSR jH 

2A+ 1 



m (IV-,A) -mo(IV+,A) 



(N(rl)-Z(r 2 p )). 

(35) 



This difference vanishes in the present case because the 
isospin impurity of the ground-state of 4 He is ignored. 



2. Energy-weighted sum rule 

Now we discuss the EWSR for the SD operator. The 
SD EWSR can be derived in the same manner as the El 
operator, and it is expressed as 



mi(p, A) = / ES(p,X,E)dE 



(36) 



where X?^ Q (H) denotes the double commutator of the 
Hamiltonian with the SD operator 



(37) 



The double commutator of the kinetic energy operator 
T = Yli=i Ti — T cm is worked out in Appendix C. The 
commutator was considered in Ref. [501 for IS and IVO 
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cases. The result for all SD cases is summarized as 



_ (A - i)h 2 



'• jv 



(2A + 1)N P 



QArriN 



(2A+1) £ {Vi-Vi)^ 



2>t=l 
A 



6m N 



6m at 



(2A + l)^(p J .(p l -iP tot ))[Tf t ,Tf] 
C^^xfp,-^)).^), (38) 



i=l 



where J-tot = Pi ls the total momentum and C p is 

related to C\ of Eq. (C7) as 



^is _ ^iivo _ r< n LV + — r«v- _ ^ r> 



(39) 



and N p = Y,t=i T ? T f reduces to A for p = IS, IVO, 
y4 — Z for p = IV+, and A — N for p = IV—, respectively. 
The isospin commutator \Tf\Tf\ vanishes for p = IS 
and IVO, while it reduces to T r z(*) for p — IV±. The 
round bracket (a x b) stands for the vector product of a 
and b, (a x 6) M = —>/2i[a x 

We name the four terms on the right-hand side of 
Eq. (|3"B"|) as model-independent (MI), spin-spin (SS), dila- 
tion (DL), and spin-orbit (SO) terms, respectively. The 
name of dilation is adopted because (p^ • (pi — ^P tot )) 
is a generator for the dilation operator. The MI term 
makes a contribution to the SD EWSR, independently of 
the ground-state wave function. Thus the kinetic energy 
contribution to the EWSR becomes model-independent 
in so far as the contribution of the other terms can be 
neglected compared to the MI term. For a fixed p the 
A-dependence of each term is simply given by 2A + 1 ex- 
cept for the SO term, which changes according to the 
ratio of 2 : 3 : (—5) for A = 0, 1, 2. On the other hand, 
for a fixed A the p-dependence of the four terms is a lit- 
tle complicated. The MI term changes in proportion to 
A : A : A - Z : A - N, while the SO term is in ratio of 
1:1: 1/2 : 1/2 for p = IS, IVO, IV+, IV-, respectively. 
The DL term identically vanishes for p = IS and IVO, 
and furthermore it turns out to have no contribution to 
the EWSR even for p = IV± because no isospin mixing 
is taken into account in our ground state of 4 He. 

Table IIVI lists the values of mi (p, A) together with 
the contributions of the kinetic energy term and its four 
terms to the EWSR calculated using the AV8'+3NF and 
G3RS+3NF potentials. The EWSR slightly depends on 
the potential models particularly for the IS SD strengths. 
Even in those cases the contribution of the kinetic energy 
to the EWSR remains almost the same. The contribution 
of the MI term to (X P X ^ Q (T)) is found to be more than 
74 % for all the cases, and really occupies a main por- 
tion of the kinetic energy contribution. The two interac- 
tions give almost the same contribution for the SS terms. 



Though the SO terms show some dependence on the in- 
teractions, the kinetic energy contributions (X?^ (T)) 
are found to be approximately model-independent. 

The enhancement of the computed sum rule (|36[) com- 
pared to (AT^ (T)) indicates the contribution of the po- 
tential energy to the EWSR. The enhancement factor 
for the El operator is 1.0 — 1.1 for the present nuclear 
forces 19]. The AV8' potential has a stronger tensor 
component than the G3RS potential. Because of this the 
tensor potential (SijTi ■ Tj) of the AV8' potential gives 
the larger contribution to the El EWSR. In the SD case, 
however, the enhancement is more complicated and de- 
pends on both multipolarity A and isospin label p. To 
elucidate this further, we have to calculate the double 
commutator for each piece of the nucleon-nucleon poten- 
tial as in the kinetic energy and evaluate its ground-state 
expectation value. 



IV. CONCLUSIONS 

We study both isovector and isoscalar spin-dipole (SD) 
strength functions in four-body calculations using realis- 
tic nuclear forces. Two different potentials are employed 
to see the sensitivity on the D-state probability produced 
by the tensor correlation. The SD excitation is built on 
the ground state of 4 He that is described accurately with 
use of explicitly correlated Gaussian bases. The contin- 
uum states including two- and three-body decay channels 
are discretized in the correlated Gaussians with aid of the 
complex scaling method. 

Experimental data that can directly be compared to 
the calculation are presently only the resonance param- 
eters of the negative-parity levels of A — 4 nuclei. Both 
the resonance energies and widths deduced from the SD 
and electric-dipole strength functions or the eigenvalues 
of the complex-scaled Hamiltonian are all in fair agree- 
ment with experiment. This success is never trivial con- 
sidering that most of the resonances among 15 levels have 
broad widths larger than 5 MeV. A combined use of both 
complex energies and appropriate strength functions pro- 
vides us with a robust tool to determine resonance pa- 
rameters. 

The non energy-weighted sum rule (NEWSR) of the 
SD strength function is investigated by relating it to the 
expectation values of three scalar products of the space- 
space and spin-spin tensors with respect to the ground 
state of 4 He. It turns out that our model space satisfies 
the NEWSR for each SD operator perfectly. The tensor 
operator of rank 2, Q%)oi * s sens itive to the D-state cor- 
relation in the ground state induced by the tensor force, 
and it is mainly responsible for distorting the ratio of the 
NEWSRs for the multipolarity A = 0,1,2 from the un- 
corrected ratio of 1 : 3 : 5. An experimental observation 
of this ratio is desirable since it may lead us to reveal the 
degree of tensor correlations in the ground state. The 
energy-weighted sum rule (EWSR) for the SD operator 
is also examined. A formula is derived to calculate the 
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TABLE IV: Energy-weighted sum rules of the SD strength functions, mi(p, A), in units of fm 2 MeV, calculated from different 
models for the nucleon-nucleon potentials. Contribution of each term of the kinetic energy to the sum rule is also listed. See 
text for the details. 

AV8'+3NF 

IS IVO IV+ IV- 

A = A = 1 A = 2 A = A = 1 A = 2 A = A = 1 A = 2 A = A = 1 A = 2 
m 1 (p, A) 126 782 949 218 450 766 110 227 389 109 225 383 

(Xf x (T)} 74.4 227 392 78.2 239 411 39.1 119 205 39.1 119 205 



MI 62.2 187 311 62.2 187 311 31.1 93.3 156 31.1 93.3 156 

SS 14.9 44.6 74.3 18.7 55.9 93.2 9.32 27.9 46.6 9.32 27.9 46.6 

DL 0.00 0.00 0.00 0.00 0.00 0.00 

SO -2.62 -3.92 6.54 -2.62 -3.92 6.54 -1.31 -1.96 3.27 -1.31 -1.96 3.27 



mi(p, A) - (A ( P A)0 (T)) 51.5 555 557 139 211 355 71.3 108 183 70.1 106 178 

G3RS+3NF 

IS IVO IV+ IV- 

A = A = 1 A = 2 A = A = 1 A = 2 A = A = 1 A = 2 A = A = 1 A = 2 

mi(p, A) 111 697 843 202 426 723 104 216 370 102 213 363 

(X? X)Q (T)} 73^2 227 403 76.3 236 419 38.1 118 209 38.1 118 209~ 



MI 62.2 187 311 62.2 187 311 31.1 93.3 156 31.1 93.3 156 

SS 16.0 47.9 79.8 19.0 57.1 95.1 9.51 28.5 47.6 9.51 28.5 47.6 

DL 0.00 0.00 0.00 0.00 0.00 0.00 

SO -4.97 -7.45 12.4 -4.97 -7.45 12.4 -2.48 -3.73 6.21 -2.48 -3.73 6.21 



mi(p,A) - (X? X)0 (T)) 37.8 470 439 126 189 304 65.5 97.8 160 63.5 94.6 153 



contribution of the kinetic energy to the EWSR. The dif- 
ference between the EWSR and the kinetic energy con- 
tribution shows some dependence on A as wells as the 
isospin character of the SD operator. Further study is 
needed to clarify the origin of its dependence by analyz- 
ing the contribution of each piece of the nuclear potential. 

Other T = resonances with Q — , 1 _ , 2~ and 1 + , 2+ 
exist in 4 He above and below the 2n + 2p threshold. It 
would be interesting to investigate these levels by the 
isoscalar SD excitation and some appropriate excitations 
produced by e.g., isoscalar quadrupole, magnetic dipolc, 
and spin-quadrupole operators with further attention be- 
ing paid to d + d type configurations. 

The SD strength functions are important inputs for 
evaluating neutrino-nucleus reaction rates. A calculation 
of neutrino- 4 He reaction rate is in progress as a conse- 
quence of the present study. It is desirable that the pre- 
dicted SD strength functions are tested with experimen- 
tal measurements in order for such reaction rate calcula- 
tion to be precise. 
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Appendix A: Multipole decomposition of the 
spin-dipole non energy-weighted sum rule 

Here we derive Eqs. (|26|) and (|2T[) by decomposing 
the operator £^ ^Au^Au mul tipoles. Substituting 
Eq. (JTJ) in &\ 0\ and recoupling the coordinate and 
spin operators, we obtain 

= (-i) A E ([« x ^ • x °iW T ? T I 

= E^2(«)o- (Al) 

The isospin operator Tf^T^ reads 

1, T z (i)r z (j), (*i • t s ) - t z (i%U) ± i(ti X tj) z (A2) 

for p — IS, IVO, and IV±, respectively. The coefficient 
Uxk is expressed by unitary Racah coefficients U as 



Ux K = (-l)y^i|t/(llll,A K ), (A3) 
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or more explicitly 




(A4) 



The operators Q p ^q and Q p ^q R have non- vanishing con- 
tributions only for n = and 2. It is easy to see that the 
2(2)0 term contains the tensor operator S%j. 



where both row and column labels, A and k, are arranged 
in order of 0, 1, and 2. The inverse of the matrix (U\ K ), 




(A5) 



is used to obtain the expectation value of Q P K ^ with re- 
spect to the ground state as discussed in Sec lIIIDTl Sec 
Eq. (US]). 

The multipole operator Q P K ^ consists of one- and two- 
body terms 



O p 



Q 



(«)o 



Q 



P(2) 
(*)0 



(A6) 



as shown in Eq. ([32]) . The two-body term with n = 2 is 
of particular interest because it contains the tensor op- 
erator. To see this, it is convenient to rewrite in 
terms of the relative and center-of-mass coordinates of 
two nucleons rather than the single-particle like coordi- 
nates, pi and pj . By introducing the coordinates ry and 
Rij by 



Ri 



Pj 



XA, 



2(1)0 * s decomposed to three terms: 



Q 



■P(2) 

(«)o 



P(2)r , Q p{2)R 



d «i y (i)o ' 



(A7) 



(A8) 



Appendix B: Calculation of the matrix elements of 
quadratic spatial tensors 



In this appendix, we give a formula of calculating the 



matrix element of Q P K ^ , Eq. (|26[) . The spin-isospin part 
can easily be evaluated in our spin and isospin func- 
tions, Eq. (ITT]) , so that we focus on the matrix ele- 
ment of the spatial part. As is clear from Eqs. ((32)) 
and (|A9|) . the spatial tensor operators have the form 
[a x b] K fj,, where a and b are vectors that represent one 
of the various coordinates, pi, pj, ry, Rij. It is useful 
to note that any of these coordinates can be expressed 
as a linear combination of the relative coordinate set x: 

a = YstLi 1 ^i x i = uix, and b = J]^ 1 & x i = ( x > where 
uj and £ are both (N — l)-dimensional column vectors. 
Therefore it is sufficient to show how we can evaluate the 
quadratic spatial tensor operators, [Uix x £x] KfJl , with the 
basis functions (fTg|). A detailed method of evaluation is 
presented in Ref . [34| , and here we follow its formulation 
and notation. 



First we calculate the matrix element between the gen- 
erating function 



where 



j>i=i 



g(s, A. x) = exp ^ — — xAx + sx^j , 



(Bl) 



n p(2)rR 
«(1)0 



3>i=l 
A 



Y, (hiX^lrkiX^i)^. (A9) 

j>i=l 



where s is an (N — l)-dimensional column vector whose 
ith element is a 3-dimensional vector s.j, and sx is a 

short-hand notation of Yli^i 1 s i ' x i ■ As given i n Ref . [13 > 
it reads 



(g(s',A',x)\ [Qx x Ca;] KM \g(s,A,x)) 



(27T) 

dctS 



AT-1 \ 2 



(B2) 
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where Tr stands for a trace and 



where 



B = A + A', v = s 



(B3) 



Using the (N — l)-dimensional column vector m specify- 
ing the GV we express s and s' as s — AieiUi + X 2 e 2 u 2 
and s' — \3e3ii3 + where a unit vector e; 

((ej • &i) = 1) and a parameter Xi are introduced to ma- 
nipulate the calculation of the sought matrix element. 
See Ref. [34j for details. The second term in the curly 
bracket and the exponential function of Eq. (|B2[) is sim- 
plified to 



C 2 



- Kfl. 

J2i<j PijK^je t -e 3 



vB 



(B4) 



Pij 



u,,B~ 



fi = CjB l m, gj = QB~ 



(B5) 



Here the arrow symbol indicates that both sides are 
equal as long as the calculation of the sought matrix 
element is concerned. That is, any terms that have 
Xf(ei ■ &i) — Xf dependence make no contribution to the 
matrix element, so that they can be dropped. 



(i) k — case 

In this case the term \CjB~ 1 v x (B~ 1 v]qq produces the 
same structure, with respect to XiXj(ei-ej), as the kinetic 
and mean square distance operators. See Appendix B.2 
of Ref. [3J| . The matrix element is 



(F(l 3 l 4 )lm(u3,u 4 , A' ,x)\ [uix x (x} 00 \F {LlL2)LM (u 1 ,u 2 ,A,x)) 

= -^={3^(5-^0 + J2(fi9j + fi9i)-£-} (F iW LM(u 3 , Ui, A', x)\F {LlL2)LM ( Ul ,u 2 , A, x)) . (B6) 



Compare this expression with Eq. (B.17) f34j . 
A formula for the overlap matrix element, 
(% 3 i 4 )iM("3)«4, A' ± x)\F (LlL2)LM (u 1 ,u 2 ,A, x)), is 
given in Eq. (B.10) [34j . 



(ii) k = 1 case 

The k = 1 case can be evaluated in exactly the same 
way as the spin-orbit matrix element of Ref. 34]. The 
result is 



(F(L 3 L 4 )L<M>(U3,U4,A',X)\ [u)X X fx]^ |-F(l 1 L 2 )Lm("1, u 2, A, x)) 
47T (-l\Li+L 2 +L+L' 



J2 V2L+lZ 2 (lL 1 L 2 L 3 L i L, LL'; kl) {F (UU)lfI (u 3 , U4, A', x)\F (Zl i 2 )Zm{ui,u 2 , A,x)) . (B7) 



r 



Compare this expression with Eq. (B.54) [34[. The Comparing this expression with Eqs. (B.41) and (B.42) 

barred angular momentum labels Li and L' follow the and using Eq. (B.49) [Hj], we obtain the matrix element 

definitions in Ref. [34|. The coefficient Z 2 is defined in as follows: 
Eq. (B. 48) [H. 

(iii) k = 2 case 

In this case we note that 



CbB^v x (B^v 



2/i 



4tt 



-r ^2(fi9j + /...'/< A, A , [e 4 x e 



(B8) 



i<3 
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(F(l 3 u)L'm>(u3,U4,A',x)\ [Cjx x C,x] 2l , \F {LlL2)LM (u 1 ,u 2 ,A,x)) 

= { ~%T + T' d^'jo vg{^gA» ( -i)^ (n^) 

x E + lZi(2LiL 2 I 3 i4i, iL'; fc) {i^x^M^, u 4 , A', x)\F {LlL2)La ( Ul ,u 2 , A, x)) 

L 

+y e (/^+/^)(-i) £i+£2 (nf^j 

l>k=l Li ) 



E V2L + lZ 2 (2L 1 L 2 L 3 L 4 L, LL'- hi) {F {L3L4)LAl (u 3 , u 4) A' 



(B9) 



The coefficient is denned in Eq. (B. 46) [34 



Appendix C: Contribution of the kinetic energy to 
the spin-dipole energy-weighted sum rule 



The aim of this appendix is to derive Eq. (|38|) . Intro- 
ducing an abbreviation 



v\^(i) = [pi x CTj] AM 



(CI) 



and T = J2i=i li — T cm , we calculate Xf^ (T) from the 
following expression 

=lEER^^)]pf.3f] 

m *>i=i 

+ [<(i),[T, WAM (i)]]i;Pi;? t }. (C2) 



Here use is made of the relation [AB, CD] = AC [5, £>] + 
[A, C]DB provided that [A,D] = and [B,C] = 0. The 
first term in the curly bracket is contributed only by i = j 
terms because [Tj^,Tf] vanishes for i ^ j. Using the 
commutation relation 



ih 

[T,v XfJ ,(i)] = [(pi - i-Ptot) x o-i]\n, 

m N 

we obtain the first term as 



(C3) 



First term = 



ih 
2m n 



E (\Pi x ^ ■ [(P* - 3 P tot) x ai] x )[Tf\Tf\ 



(C4) 



The ground-state expectation value of this term is conve- 
niently evaluated by decomposing the above scalar prod- 
uct to that of the space-space and spin-spin terms using 
the matrix U of Eq. (|A4|l . The result is 

First term = £ ( Pi • (p, - £p tot )) [Tf t,T?] 



x (p* " i p *ot)) ■ ^1^,^]. (C5) 



The matrix element of the spatial part involving the op- 
erators, (pi • (pi - ^Ptot)) and {(pi x (p s - ^Ptot)) • Cj)' 
can be calculated in the manner similar to that presented 
in Appendix B. See Ref. [13] for the details. 

The second term in the curly bracket of Eq. (|C2[) can 
be obtained in a similar way. After a straightforward 
calculation of the commutator, we obtain the following 
result: 



Second term 

h 2 



2th 



(2A + l)N p 



QArriN 



(2A + l)(£ p -£ pt ) 



6m 



N 



Cx E ((* x (P* - i P '°t)) • ffi)?? Tf f , (C6) 



where Ca is 

Co = 2, Ci=3, Cfe = -5. (C7) 
Here the operators N p and S p are defined by 



(C8) 



i=i 



which leads to (£p • = 3A" 5 + £f >i=1 (o-j ■ cr j)^- 
Combining (|C5[) and (|(J6|) we obtain Eq. 
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